Constraining dark energy with SNe la and large-scale structure 
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Measurements of the distances to SNe la have produced strong evidence that the expansion of the 
Universe is accelerating, implying the existence of a nearly uniform component of dark energy with 
negative pressure. We show that constraints to this mysterious component based upon large-scale 
^\ I structure nicely complement the SN la data, and that together they require Q.x £ (0.6, 0.7) and 

^\ . wx < —0.6 (95% cl), for the favored flat Universe. Other cosmological data support this conclusion. 

0^ ' The simplest explanation, a cosmological constant, is consistent with this, while some of the other 

possibilities are not. 
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I. INTRODUCTION 



Two groups [gJJ^ have presented strong evidence that the expansion of the Universe is speeding up, rather than 
slowing down. It comes in the form of distance measurements to some fifty supernovae of type la (SNe la), with 
, redshifts between and 1. The results are fully consistent with the existence of a cosmological constant (vacuum 
energy) whose contribution to the energy density is around 70% of the critical density (JIa ^ 0.7). Other measurements 
indicate that matter alone contributes VLm = 0.4 ±0.1 |^]. Taken together, matter and vacuum energy account for 
an amount close to the critical density, consistent with measurements of the anisotropy of the cosmic microwave 
background (CMB) g. 

' In spite of the apparent success of the cosmological constant explanation, other possibilities have been suggested 
for the "dark energy." This is in part because of the checkered history of the cosmological constant: It was advocated 
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O ■ by Einstein to construct a static universe and discarded after the discovery of the expansion; it was revived by Hoyle 
^ and Bondi and Gold to solve an age crisis, later resolved by a smaller Hubble constant, and it was put forth to explain 
^ the abundance of quasars at 2 ~ 2, now known to be due to galactic evolution. Further, all attempts to compute the 
. . , value of the cosmological constant, which in modern terms corresponds to the energy associated with the quantum 
. ^ vacuum, have been wildly unsuccessful ||^ . Finally, the presence of a cosmological constant makes the present epoch 
, special: at earlier times matter (or radiation) dominated the energy density and at later times vacuum energy will 
' dominate (the "why now?" problem). 

The key features of an alternative form for the dark energy are: bulk pressure that is significantly negative, 
w < —1/3, where w = p/p, and the inability to clump effectively. The first property is needed to ensure accelerated 
expansion and to avoid interfering with a long matter-dominated era during which structure forms; the second 
property is needed so that the dark energy escapes detection in gravitationally bound systems such as clusters of 
galaxies. Candidates for the dark energy include |^: a frustrated network of topological defects (such as strings or 
walls), here w = (n is the dimension of the defect) |^] and an evolving scalar field, where p — -f V{(t>) and 

p = — V{(j)) (referred to by some as quintessence) |^,^. 

The SN la data alone do not yet discriminate well against these different possibilities ||l|,|l^. As shown in Fig. 
the maximum likelihood region in the ^m~w plane runs roughly diagonally: less negative pressure is permitted if 
the fraction of critical density contributed by dark energy is larger. Following earlier work this led us to consider 
other cosmological constraints: large-scale structure, anisotropy of the CMB, the age of the Universe, gravitational 
lensing, and measurements of the Hubble constant and of the matter density. As we shall show, some of the additional 
constraints, especially large-scale structure, complement the SN la constraint, and serve to sharpen the limits to ^Im 
and w; others primarily illustrate the consistency of these measurements with the SN la result. In the end, we find 
nx e (0.6,0.7) and w < -0.6 (95% cl). 
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II. METHOD 



Our underlying cosmological paradigm is a flat, cold dark matter model with a dark-energy component, though as 
we will discuss later our results are more general. We restrict ourselves to flat models both because they are preferred 
by the CMB anisotropy data and a flat Universe is strongly favored by inflation. We restrict ourselves to cold dark 
matter models because of the success of the cold dark matter paradigm and the lack of a viable alternative. For our 
space of models we construct marginalized likelihood functions based upon SNe la, large-scale structure, and other 
cosmological measurements, as described below. 

Our model parameter space includes the usual cosmological parameters (f^M, ^sh^, and h) and the amplitude 
and spectral index of the spectrum of Gaussian curvature fluctuations (cts and n). For the dark-energy component, 
we choose to focus on the dynamical scalar-field models, because the frustrated defect models are at best marginally 
consistent with the SN la data alone . 

In the dynamical scalar-field models the equation of state w = p/ p varies with time. However for most of our 
purposes, only one additional free parameter needs to be specified, an "effective" equation of state. We choose Wcb 
to be that value w which, if the Universe had w constant, would reproduce the conformal age today. We choose 
this definition because the CMB anisotropy spectrum and the COBE normalization of the matter power spectrum 
remain constant (to within 5-10%) for different scalar field models with the same Woff 0|. For the models under 
consideration vues is closely approximated by 

Wcff = J da V.^{a)w{a)/ J da Vt^{a) . (1) 

and, since it is simpler to compute, we have used Wcs throughout. Obviously, our results also apply to constant w 
models (e.g., frustrated defects), by taking w — WcG- 

While WeS neatly parameterizes the scalar-field models from the standpoint of large-scale structure and the CMB 
anisotropy, it does not do as well when it comes to the SN la data. Recall that w^s as defined in Eq. (|l|) receives 
a contribution from a wide range of redshifts. The SN la data however are sensitive mostly to z ^ 1/2. Since w 
becomes less negative with time in the models we are considering, the SN la data "see" a less negative w than the 
CMB by a model dependent amount. We shall return to this point later. 

We normalize our models to the COBE 4-year data |l^ using the method of Ref . |l^ . Beyond the COBE measure- 
ments, the small-scale anisotropy of the microwave background tells us that the Universe is close to being spatially 
flat (position of the first acoustic peak) and that VLm is less than one and/or the baryon density is high (height of 
the first acoustic peak). We have not included a detailed fit to the current data (see e.g. Ref. Q), but rather impose 
flatness. The additional facts that might be gleaned from present CMB measurements, J^m < 1 and high baryon 
density, are in fact much more strongly imposed by the large-scale structure data and the Buries - Tytler deuterium 
measurement. 

We require that the power-spectrum shape fit the redshift-survey data as compiled by Ref. [p^ (excluding the 
4 smallest scale points which are most sensitive to the effects of bias and nonlinear effects). On smaller scales we 
require that all of our models reproduce the observed abundance of rich clusters of galaxies. This is accomplished by 
requirmg dg = (0.55 ± 0.1) where CTg is the rms mass fluctuation in spheres of 8 /i ^ Mpc computed in linear 

theory [p^ . The baryon density is fixed at the central value indicated by the Burles-Tytler deuterium measurements, 
ilsh^ = 0.019 ±0.001 ||l^. We assume that clusters are a fair sample of the matter in the Universe so that the cluster 
baryon fraction Jb = (0.07 ± 0.007)ft.~^/^ reflects the universal ratio of baryons to matter {Qb/^m)- We marginalize 
over the spectral index and Hubble constant, assuming Gaussian priors with n — 0.95 ±0.05, which encompasses most 
inflationary models, and h = 0.65 ± 0.05, which is consistent with current measurements. 

There are three other cosmological constraints that we did not impose: the age of the Universe, — (14 ± 2)Gyr 
[ p^ ; direct measurements of the matter density, Om = 0.4 ±0.1, and the frequency of multiply imaged quasars. 
While important, these constraints serve to prove consistency, rather than to provide complementary information. 
For example, the SN la data together with our Hubble constant constraint lead to an almost identical age constraint 
[0,^. The lensing constraint, recently studied in detail for dynamical scalar-field models |l9j, excludes the region of 
large fix and very negative w (at 95% cl, below the line WcS = —0.55 — I.S^Im), which is disfavored by the SN la 
data. The matter density determined by direct measurements, ^Im = 0.4 ±0.1, is consistent with that imposed by 
the LSS and Hubble constant constraints. 
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III. RESULTS 




FIG. 1. Contours of likelihood, from 0.5a to 2a, in the Qm-WeS plane. Left: The thin solid lines are the constraints from 
LSS and the CMB. The heavy lines are the SN la constraints (using the Fit C supernovae of Ref. for constant w models 
(solid curves) and for a scalar-field model with an exponential potential (broken curves; quadratic and quartic potentials have 
very similar SN la constraints). Note that the SN la contours for dynamical scalar-field models and constant w models are 
slightly offset (see text). Right; The likelihood contours from all of our cosmological constraints for constant w models (solid) 
and dynamical scalar-field models (broken). 

As can be seen in Fig. |l|, our large-scale structure and CMB constraints neatly complement the SN la data. LSS 
tightly constrains ^m, but is less restrictive along the Wcs axis. This is easy to understand: in order to fit the power 
spectrum data, a COBE-normalized CDM model must have "shape parameter" F — Vt^h ^ 0.25 (with a slight 
dependence on n). Together with the constraint h = 0.65 ± 0.05 (and our Jb constraint) this leads to Hm ~ 0.35. 
As discussed in Rcf. the erg constraint can discriminate against Wcs', however, allowing the spectral index to differ 
significantly from unity diminishes its power to do so. 

Note that the SN la likelihood contours for the dynamical scalar-field model and the constant-w models are not 
the same while the LSS contours are identical. With the Fit C supernovae of Ref. 0] and the dynamical scalar-field 
models considered here (quadratic, quartic and exponential scalar potentials), the contours are displaced by about 
0.1 in WcS'. the 95% cl upper limit to WcS for the constant w models is —0.62, while for the quartic, quadratic and 
exponential potentials for V{(l)) it is —0.75, —0.76 and —0.73 respectively. The reason for this shift is simple: the w 
dependence of LSS is almost completely contained in the distance to the last-scattering surface and Wcs is constructed 
to hold that constant. On the other hand, the w dependence of the SN la results is more heavily weighted by the 
recent value of w; said another way, there is a different effective w for the SN la data. This fact could ultimately prove 
to be very important in discriminating between different models. 

Additionally there are a class of dynamical scalar-field models that have attracted much interest recently [p|,pO|. 
For these potentials (here we consider V{(j)) — c/^P and V{(l)) — c[e^/"^ — 1]), and a wide range of initial conditions 
the scalar-field settles into a "tracking solution" that depends only upon one parameter (here c) and the evolution of 
the cosmic scale factor, suggesting that they might help to address the "why now?" problem. 

For our purposes, the most interesting fact is that each tracker potential picks out a curve in ^y/ — ^cff space. 
Typically the lower values of VLm go with the most negative values of Wcs and vice versa (see Fig. H) This fact puts 
the tracker solutions in jeopardy, as shown in the same figure. For the tracker models shown here (p = 2, 4 and 
exponential), the 95% cl intervals for the SN la and LSS data barely overlap. The situation is even worse for larger 
values of p. A similar problem was noted in Ref. 

Finally, we comment on the robustness of our results. While we have restricted ourselves to flat models, as preferred 
by the CMB data, our constraints do not depend strongly on this assumption. This is because the LSS constraints 
are insensitive to the flatness assumption, and curvature, which corresponds to a Woff = —\ component, is strongly 
disfavored by the SN la results. We have not explicitly allowed for the possibility that inflation-produced gravity 
waves account for a significant part of the CMB anisotropy on large-angular scales (i.e., T/S > 0.1), which would 

3 



have the effect of decreasing the overall amplitude of the COBE normalized power spectrum. In fact, allowing for 
gravity waves would not change our results, as this degree of freedom is implicitly accounted for by a combination of 
n, the normalization freedom in the power spectrum and the uncertainty in the COBE normalization. 

Our model space does not explore more radical possibilities, for example, that neutrinos contribute significantly to 
the mass density or a nonpower-law or isocurvature spectrum of density perturbations p^ . Even allowing for these 
possibilities (or others) would not change our results significantly if one still adopted the mass density constraint, 
~ 0.4 ± 0.1. As discussed earlier, it is almost as powerful as the CDM-based LSS constraint. 
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FIG. 2. Upper panel: The relationship between WeS and Qm for a selection of tracker potentials. Lower panels: the CMB 
and LSS likelihoods from Fig. ^ as a function of (dotted) and the SN la likelihood (solid - normalized to unity at the 
peak). As can be seen clearly, tracker models have difficulty simultaneously accommodating the SN la and LSS constraints. 



IV. CONCLUSIONS 



The evidence provided by SNe la that the Universe is accelerating rather than slowing solves one mystery - the 
discrepancy between direct measurements of the matter density and measurements of the spatial curvature based 
upon CMB anisotropy - and introduces another - the nature of the dark energy that together with matter accounts 
for the critical density. SNe la alone do not yet strongly constrain the nature of the dark energy. 

In this Letter we have shown that consideration of other important cosmological data both complement and reinforce 
the SN la results. In particular, as illustrated in Fig. [l], consideration of large-scale structure leads to a constraint 
that nicely complements the SN la constraint and strengthens the conclusions that one can draw. Other cosmological 
constraints - age of the Universe, frequency of gravitational lensing and direct measures of the matter density - 
provide information that is consistent with the SN la constraint (lensing and age) and the LSS constraint (matter 
density) , and thereby reinforces the self consistency of the whole picture of a flat Universe with cold dark matter and 
dark energy. 

Finally, what have we learned about the properties of the dark-energy component? The suite of cosmological 
constraints that we have applied indicate that fix G (0.6, 0.7) and WcS < --0.6 (95% cl), with the most likely value of 
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WeS close to —1 (see Fig. p. The frustrated network of light cosmic string {wcs = —5) is strongly disfavored, and a 
network of frustrated walls {wcs = — |) is only slightly more acceptable. Also in the disfavored category are tracker 
models with V{(j)) = c/(ff and p = 2, 4, 6, 8, • • •. Dynamical scalar-field models can be made acceptable provided 
is tuned to be more negative than —0.7. The current data definitely prefer the most economical, if not the most 
perplexing, solution: Einstein's cosmological constant. 
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